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We study a simple system described by a 2 x 2 Hamiltonian and the evolution of the quantum 
states under the influence of a perturbation. More precisely, when the initial Hamiltonian is not 
degenerate, we check analytically the validity of the adiabatic approximation and verify that, even if 
the evolution operator has no limit for adiabatic switchings, the Gell-Mann and Low formula allows 
to follow the evolution of eigenstates. In the degenerate case, for generic initial eigenstates, the 
adiabatic approximation (obtained by two different limiting procedures) is either useless or wrong, 
and the Gell-Mann and Low formula does not hold. We show how to select initial states in order to 
avoid such failures. 
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I. INTRODUCTION 

Adiabatic switching is a crucial ingredient of many- 
body theory. It provides a way to express the eigenstates 
of a Hamiltonian H = Hq+Hi in terms of the eigenstates 
of Hq. Its basic idea is to switch very slowly the inter- 
action Hi, i.e. to transform H into a time-dependent 
Hamiltonian H{t) = Hq + c^^'^^'Hi where the small pa- 
rameter £ > eventually vanishes. Under the influence of 
H{t), an eigenstate |To) oi Hq becomes a time-dependent 
wavefunction [^^(i)) and it might be expected that an 
eigenstate oi H — H{0) is obtained by taking the limit 
of l^e(O)) when e tends to zero. It turns out that this 
naive expectation is not justified because |\I'£(0)) has no 
limit when e ^ 0. When the initial state |To) belongs 
to a non degenerate eigenspace (isolated from the other 
eigensubspaces), Gell-Mann and Low [1| solved the prob- 
lem by dividing out the divergence by a suitable factor. 
The ratio is called the Gell-Mann and Low wavefunction 
and its convergence can be proved by using the adiabatic 
theorem [2]. 

In the first part of this work, we present an exactly 
solvable 2x2 model that illustrates the fact that the limit 
e ^ of |*£(0)) does not exist. The validity of the Gell- 
Mann and Low wavefunction is shown by analytically 
calculating the corresponding adiabatic approximation. 

It was realized fifty years ago [3| that the Gell-Mann 
and Low formula must sometimes be extended to the 
case of a degenerate initial state of Hq. This happens in 
many practical situations, for instance when the system 
contains unfilled shells. This problem has been discussed 
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in several fields, including nuclear physics [3], solid state 
physics , quantum chemistry Q and atomic physics [3] . 
In most cases, it is assumed that there is some eigenstate 
I To) for which the Gell-Mann and Low formula holds. 
More precisely, if Vo is the vector space generated by the 
eigenstates of Hq associated with the degenerate ground 
state energy Eq, then the claim is that there exists some 
initial state |T) in Vo whose time evolved state |^e(0)) 
is (up to a divergent phase) an eigenstate of H. To see 
that it is not possible to choose any element of Vo as ini- 
tial state, let us forget for a moment the divergent phase 
and consider a perturbation Hi that splits the degen- 
eracy and two initial states |T^) and |T^) whose time 
evolutions give rise to two eigenstates of H, denoted by 
1*^(0)) and 1*^(0)), with different energies Ei and £'2. 
By linearity, the initial state jT-'^) -f- |T^) evolves into 
1*^(0)) -I- 1*^(0)) which is not an eigenstate of H be- 
cause El ^ E2. Therefore, |T^) + |T^) is not a proper 
initial state. 

When the initial state of Hq is degenerate, we show 
that, even for the very simple 2x2 model considered 
here, the Gell-Mann and Low wavefunction does not have 
a limit for almost all initial conditions. We find however 
the specific initial states that lead to convergent wave- 
functions. For the application of many-body methods to 
degenerate systems, it is crucial to find a way to select 
the proper initial states. 



II. NON-DEGENERATE INITIAL STATES 

Two-dimensional matrices are the simplest non-trivial 
models that can be considered in quantum physics, and 
were indeed already used as toy-models for many-body 
theory studies (see ReL [8| and references therein) . How- 
ever, it is the first time to our knowledge that the evolu- 
tion operator and its adiabatic approximation are calcu- 



lated explicitly for such a model. 

We consider the two-dimensional system described by 
the Hamiltonian H = Hq + Hi with 



Hn = 



^1-5 
fi + S 



Hi = 



X 
X 



(1) 



We assume that (5 7^ 0, so that the initial state is non de- 
generate. Without restriction, it can further be assumed 
that 6>Q. 

According to the usual treatment of adiabatic switch- 
ing [9|, llO| , it is convenient to transform the Schrodinger 
equation for the time-dependent Hamiltonian H{t) — 
Hq + e^'^l*l_ffi to the interaction picture. The depen- 
dence of the operators and wavefunctions on the param- 
eter e > will be implicit in the sequel. If we denote by 
'^ s{t) a solution of the Schrodinger equation 

the wavefunction in the interaction picture is defined as 

*(t) = e'^°**s(i)- 
It satisfies 

i^-^^Hum{ti 

subject to the boundary condition $(—00) = T, where 
T may be an eigenstate of Hq. The Hamiltonian Hi^tit) 
is the Hamiltonian H{t) in the interaction picture 

Hint{t) = e-^l*l e'^°* Hi c-''"°K 
In the simple case ([T]) considered here, this operator reads 

iJint(i)=a:e-^l*l 



/ n „~2\5t 
-,-e\t\ I u e 

^2i5t Q 



(2) 



Instead of using wavefunctions, it is customary to work 
with the evolution matrix U{t) such that ^(i) = U{t)T. 
We refer to JJJ for sufficient conditions on the time de- 
pendent Hamiltonian i?int(0 i'^ order to ensure the exis- 
tence of the unitary propagator U{t). 



A. The evolution operator: analytic solution and 
limiting behavior 

We give in this section the analytic expression of the 
evolution matrix. This operator is the solution of the 
Schrodinger equation in the interaction picture 



i^ = i/,„,(i)[/(i), 



with the boundary condition U{—oo) = Id. In the se- 
quel, we consider i < 0, so that the switching function 



is e"^*. We also assume that x > 0. Denoting the matrix 
elements of U{t) by 



_ / a{t) b{t) \ 



(4) 



eq. ([31) is equivalent to the system of equations 

ia' = a;e-'(25+i£)tc, 

ic' = a;e'(2'5-i'^)*a^ 

ib' = xe-'^^^+'^'>*d, 

id' = xe'<^^^-'^'>*b, 

with the boundary conditions 

a(— 00) = d{—oo) = 1, 6(— cxd) = c(— 00) = 0. 

Since the functions c, d are independent of the functions 
a, b, and satisfy equations of the same form as for a, b, 
it is enough to solve the first two equations of the above 
system. 

1. Solution of the equation for the unknown function a 

By eliminating the function c in eq. ([?]), the second 
order evolution equation 



a + (2i5 - e)a + x^e^^^a = 



(5) 



is obtained. This equation can be solved by using stan- 
dard techniques [l2, Eq. (23), p. 442]: If we rewrite the 
unknown function a as a{t) = e~'*^''+"^/^^*Ziy(a;e'^*/£); ^'^^ 
introduce the variable s = xe^'/e, eq. (O becomes the 
Bessel equation for Z,j{s), with 



v = l/2-iS/e. 



(6) 



Since Ji, and J_j/ are two independent solutions of the 
Bessel equation when v is not an integer [iJl , the function 
a{t) = d{s) has the general form 



d{s) = (^^Y {CiMs) -f C2J-AS)) 



(7) 



where the constants Ci , C2 are determined by the bound- 
ary conditions at s = (i.e in the limit t — > —00). Since 
a(— 00) = a(0) = 1, and by the series expansion |13l . 
Eq. (9.1.10)] 



'W^(I) 



sYys (-sV4)^ 

/c=0 



it follows Ci = and C2 = 2-'"r(l - i/). 



(3) 2. Solution of the equation for the other unknown functions 



The function c{t) — c{s) can be obtained from the 
expression of a(t) — d{s), relying on the first equation 



in the system ((4]), rewritten in the s variable as c(s ) = 
i{es/xy-'^''a'{s). Since (s^J^^is))' = -s'^Ji_^(s) [H, 
Eq. (9.1.30)], it holds 



c{s) 



'IC2 



l~2iy 



(-) s^-^J.-.is). 



The functions 6 and d satisfy the same equations as a and 
c, respectively, but with different boundary conditions. 
With the above notation, it can easily be checked that h 
is also a combination of Bessel functions as given by ([7]), 
but with Ci = -i[e/xf''-^2''^^T{v) and C2 = 0. 



3. Properties of the solution 

In view of the above results, the problem ([4| has the 
analytic solution 



b{t) = Ci (^)%-*J, (" 



xe' 

e 



c{t) 



-\C2 ( - I e 



r(l-i.)t 



Jl- 



qE* 






d(i) = iCi (^)%-(i-)*j._i(^^). 



y^. Limiting behavior of the evolution operator 

Figure [1] illustrates the divergence of a(0) as a func- 
tion of e. The above analytic expressions therefore show 
that the evolution operator does not have a limit when 
e — > 0. Rather, as illustrated by Fig. [H a strongly os- 
cillatory behavior is observed. This is the reason why 
the phase factor has to be cancelled out by considering 
a renormalized wavefunction in the Gell-Mann and Low 
fashion, as explained in the next section. 
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FIG. 1: Real part of a(0) as a function of e for x — 5 — 1 and 
fj, = 0. For clarity, the region x < 0.005 is not plotted. 



with 



and where i^ is defined in ([5]). As a consistency check, it is 
possible to verify that the matrix U is unitary. This fol- 



lows from \Ci 



ICar = r(i^)F(l - i/)/2 = 7r/(2sinzy7r) 



and J^(s)Ji_i/(s) -I- J^^{s)Ji,^i{s) = 2sin(z/7r)/(7rs) (see 
[H Eq. (9.1.15)]). 
The series expansion of the Bessel functions gives 



00 



a{s) 
bis) 
cis) 
dis) 



1 + E 



i-s^/Af 






2Sfc!n-;i(j+--i)' 



JliSt ■ 



(-sV4)* 



2 S^! mil (.?--)' 



1+E 



3 = 
2 /A\k 



(-sV4) 



t[klUU(-^+'^-^) 



Note that these solutions are valid for all values oix, 5 ^ 
and e > 0. We have the symmetries a{t, —S) = d{t,S) 
and c{t,~5) = -b{t,S), a{t)* = d{t) and b{t)* = c{t). 
Moreover, a{t) and d{t) are even in x, while b{t) and c{t) 
are odd in x. 



B. The Gell-Mann and Low w^avefunction 

Rigorous proofs of the Gell-Mann and Low formula 
rely on the adiabatic approximation. Therefore, we first 
compute the adiatic approximation for the problem ([1]), 
and then the Gell-Mann and Low wavcfunctions. 



1. The adiabatic approximation 

The adiabatic approximation consists in approximat- 
ing the evolution operator in the interaction picture U(t) 
by the adiabatic operator in the interaction picture: 



with 



Uait)=e''""AiT)<^,{T) 



$,(t)= lim $,(T,To)e-"°^°, 

To— ►-00 



(8) 



where t = et < and tq = eio [iJ, \lM ■ The unitary ma- 
trices A and <&£ are calculated from the eigenprojectors 
of H{t) ~ Hq + e'^Hi, as explained below (see formulas 
(Uni) and (HID). 

The eigenvalues of H{t) are 



e± 



(r) = M ± A, A = VP + 



2p2T 



x^e 



With the notation, y = xe^ /5, the corresponding eigen- 
projectors read respectively 



P±ir) 



1 



2yTT7 



VT+Fti ±y 



±y yiTFii 



Defining the Herniitian matrix [If 



ly 



dr 



1 



_± ., 2(l + y2) V-1 Oy 

the matrix A appearing in ([8]) satisfies the equation 

A'{t) = -ii^(T)A(r), (9) 

with the boundary condition A{—oo) ~ Id. The unique 
solution is 



A(t) 



1 / vT+a \/l — a 



V2 V -\/l-a Vl + a/ ' 



(10) 



with a = (l + y2)-i/2 £ [0^1]. 

The phase matrix $e(T, tq) is obtained from the inte- 
gral of the eigenvalues: 



<I'e(r,ro) = Y. ^^P {'- r^Arldr') P,( 



cr = ± 



[ro) 



For large negative tq, up to 0(e'^'') terms coming from 
the approximation P(j(to) ~ Pa-{—oo), 



$^.(r,ro)=e~'''(*-*«) 
where 

with 



gi0(r,To)/e Q 

Q-i<t>iT.To)/e 



, (11) 



f y^sm!^da = F{T)-F{TQ), 



Note that ^^(to) has no finite limit for tq -^ — oo. How- 
ever, because the computations are done in the interac- 
tion picture, the phase matrix ^i;{t,to) is multiplied by 
the unitary operator e"''"^", and the so-obtained opera- 
tor has a limit when tq — > — cxd. Indeed, 



$£(r) = lim $£(T,ro)e"'^''*" 
_ p-i^r/e f exp(i</>(r)/e) 







exp(— i0(T)/e) 



(r) = lim ^To + (pir, tq) - (// - ^)to 

Co — > — OO 

25 



where 



= F(t) - 5 I 1 - lo: 

At i = 0, we obtain the adiabatic approximation 

.. ,„. _ If e''^(o)/^VTT^ e-''^(o)/^yr^^ 
"^ ^ " V2 V ~e'*(o)/-yT^7^ e-W)/^x/rT^ 

where a is now (1 -I- x^ /5^)~^^^. 



2. Quality of the adiabatic approximation 

The rigorous adiabatic theorem (see Ref. 16 for a recent 
account) states that the adiabatic evolution operator Ua 
defined by Eq. ([5]) is such that 

snp\\U{t)-Uam <Ce, 

for some constant C. In particular, for i = 0, ||C/(0) — 
t^a(0)|| < Ce. Figure 2 shows that this result is indeed 
verified in the case considered here, and allows to give a 
numerical estimate of the constant C. 
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FIG. 2: Modulus of the difference between a(0) and its adi- 
abatic approximation as a function of e for x = 5 — 1 and 
Ai = 0. 



3. The Gell-Mann and Low limit 

The Gell-Mann and Low formula describes how an 
eigenvector of the reference (unperturbed) Hamiltonian 
Hq evolves under an added perturbation Hi. Starting 
from one of the initial eigenstates of iJo, namely 



|Ti) = 



|T2) = 



(12) 



the associated Gell-Mann and Low wavefunctions 

C/(0)|T') ^. C/a(0)|T») 



|^X0)> = lim 



lim 



-0 (T*|C/(0)|T^) e^O (T*|C/a(0)|T*) 
can be computed using ^. They read respectively 

1 



and 






X 




VaF 


2+^2 


+ 5 


Va;2 


+ 5^ + 5 



X 

1 



where aij{T) are the matrix elements of A{t) given by 
PH)) . It is easy to check that these vectors are indeed 
ei genstate s of H, for the eigenvalues /x— ^/x'^ + 5'^ and fi+ 
\J x^ + (5^ respectively. These eigenvalues can be obtained 
by the energy-shift formula [lO, P- 200]. The difference 
Ai5 between the energy of an eigenstate of i^o + ^i and 
of the corresponding eigenstate |To) oi Ho is 



where U(t) is given by Eq. ([13]). The adiabatic approx- 
imation is therefore exact but, as we shall see, the Gell- 
Mann and Low wavefunction has no limit for the initial 
states T^ and T^ when e ^ 0. 

Taking first the limit tq — > — oo amounts to take the 
limit (5 ^ in the definition of the adiabatic evolution 
operator ([5]). For the operator A, the limit is 



A^ = limiea;— log(To|C/(0)|To) 
e^o ax 



limiea:— log(To|[/a(0)|To> 
E^o dx 



We find indeed that AE — S - 
AE=-S + Vx^ + (52 for \T^). 



Vx^ + P for |Ti) and 



III. DEGENERATE CASE 

Initial degenerate states for the model ([T]) are obtained 
when 6 — 0. In this case, the evolution operator can 
still be computed analytically, and actually has a simpler 
expression than the one obtained in Section III Al in the 
non-degnerate case: 



Uit) 



_ f cos(xe^ /e) — isin(a;e^7£) 
-isin(xe^*/£) cos(a;e^*/£) 



(13) 



•^'''^7l( 



1 1 
-1 1 



(14) 



for any r < 0. Notice that the boundary condition at r = 
— oo for the equation ([5]) on A is therefore not satisfied. 
For the matrix ^^(r), the limit is 



*e(r) 



— p-iM* 



exp(ixe'^/e) 

exp(— ixe'^/e) 



The limit of the adiabatic evolution operator reads now 

. ._ 1 f exp(ia:e^Ve) exp(-ia;e^Ve) \ , . 
^^' ~ ^/2 y-eM^xe^ys) exp(-ia;e^7e) J ' ^^^^ 

It is however not an approximation of U{t). 

The discrepancy between the two approaches clearly 
shows that the two limits 5 — > and tq — > — oo do not 
commute for degenerate systems. 



In particular. 



TTfn.\ _ ( cos(a;/e) — isin(x/e) 
^ ' \ — isin(a;/e) cos(a;/e). 



A. Failure of the adiabatic approximation and the 
Gell-Mann and Low formula 

1. Adiabatic approximation 

There are two ways to calculate the adiabatic approx- 
imation, depending on the order of the limits 5^0 and 

To -> -OO. 

If we first carry out the limit (5^0, the computa- 
tion of Section IIIB II can be repeated by starting from 
the Hamiltonian Hq with (5 = 0. The eigenvalues corre- 
sponding to H{t) are /i ± xeT , and the eigenprojectors 
are constant: 



P.{r)^\ 



1 ±1 
±1 1 



Thus, K{t) — 0, so that A{t) is constant and equal to its 
boundary value A{t) — Id. The computation of <I>e(r) is 
therefore straighforward and leads to the adiabatic ap- 
proximation 

Ua{t) = U{t), 



2. Non-validity of the general Gell-Mann and Low formula 

Consider a given initial state |To), which is an eigen- 
vector associated with the eigenvalue ^. Since the corre- 
sponding eigenspace is two-dimensional, there seems to 
be some arbitrariness in the choice of the initial state. 
However, we shall see that most initial states lead to 
divergent Gell-Mann and Low wavefunctions. If for in- 
stance |To) = |T^) = I „ j as in Section III B 31 the Gell- 
Mann and Low wavefunction 



t/(0)|To) 
(To|[/(0)|To) 



-i tan(a;/e) 



has no limit when e ^ 0. This shows that the Gell-Mann 
and Low formula fails for the initial state |T-'^). It fails 
also for |T2). 

Similarly, the energy shift would be given by the limit 
for e ^ of 

AE = \ex-^\og{T^\U{Q)\T^) = -ixta.n{x/e). 
ax 

But this has no limit when a; 7^ 0. 

The question is then whether there are some initial 
states for which the Gell-Mann and Low wavefunction 
has a limit, and how those states can be characterized. 



B. Selection of the proper initial states 

In this section, we try to find initial states that are 
eigenstates of Hq and that lead to convergent Gell-Mann 
and Low wavefunctions. 

Any initial state |To) for the model H]) in the degen- 
erate case is of the general form (up to a trivial scaling 
changing this state into — |To)) 



|Tn)=cos0|Ti) 



sm f 



where < 6* < tt and \T^),\T^) are defined in p2)) . 
Straightforward computations show that 

C/(0)|To) 



(To|t/(0)|To)' 

has a limit if and only if 6* = 7r/4 or 9 ~ in/ A. This 
defines two proper initial states 

with associated Gell-Mann and Low wavefunctions 

f/(0)|T±) 



(T±|[/(0)|Tj 



1 

±1 



which obviously have limits in the regime e = 0. More- 
over, in the basis (|T_), |T-|_)), the evolution operator 
becomes 

fjC^N ^_}_( exp(ixe^7e) exp(-ia;e^7£) \ ,^g. 
^ ' y/2 V ~ exp(ia;e'^*/£) exp(— ia;e'^*/£) ) 



Similarly, the energy-shift formula gives now the correct 
result A_E± = ±x. 

Interestingly, the adiabatic operator given in eq. (|15p is 
equal to the evolution operator in the basis (|T_), |T+)). 



IV. CONCLUSION 



The adiabatic theorem, which was recently ques- 
tioned [l3|, was investigated here by using an exactly 
solvable model. Within this model, the adiabatic theo- 
rem is valid when the initial state of the system is non- 
degenerate. When it is degenerate, our simple model 
exhibits several problems that are generally present for 
an arbitrary initial state: the Gell-Mann and Low wave- 
function does not converge, the energy-shift formula is 
not valid and the adiabatic approximation becomes am- 
biguous. Within our model, all these problems are solved 
by properly choosing the initial state. 

At least since Tolmachev [1^, it is conjectured that 
the Gell-Mann and Low formula is valid for properly cho- 
sen initial states when the system is degenerate. How- 
ever, this conjecture has not been proved and no practical 
method was given to select the proper initial states. We 
intend to come back to this question in a forthcoming 
paper. 
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